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ABSTRACT 

Using measured tangential shear profiles and number counts of massive elliptical galaxies, the halo occupation 
distribution of strong lensing galaxies is constrained. The resulting HOD is then used to populate an N-body 
simulation with lens galaxies, in order to assess the importance of environment for strong lensing systems. 
Typical estimated values for the convergence and shear produced by nearby correlated matter are K^i-p — Ygrp — 
0.03, with much stronger events occurring relatively infrequently. This implies that estimates of quantities like 
the Hubble constant are not expected to be significantly biased by environmental effects. One puzzle is that 
predicted values for the external shear at lens galaxies are far below the values obtained by modeling of strong 
lensing data. 

Subject headings: Gravitational lensing 
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1. INTRODUCTION 

Strong gravitational lensing has emerged as an impor- 
tant tool in the astrophysical toolbelt. Strong lenses have 
been u sed, for example, to measure the Hubble constant (e.g 
ISchechter..2004 . and references therein), to study the proper- 
ties of elliptical galaxies (e.g. Rusin et al. 2003), to study dark 
matter substructure (e.g. Kochanek & Dalal 2004), to resolve 
the inner structure of quasars (Richards et al. 2004), even 
to co nstrain the geometry of the universe (Broadhur st et at] 
120041). Many of the applications of strong lensing require re- 
constructing the gravitational potential of the lens, which in- 
cludes contributions both from the principal lens galaxy and 
from any (projected) nearby structures. Models of strong 
lenses often treat the lens galaxy itself in exquisite detail (e.g 
ICohn et al. 2001) . however the lens environment is usually 
treated more simply. Except in ca s es with obvious n earby 
galaxies (e.g.' Koopmans et alJl2()03t iKundic et alJll997 ). the 
environments of lenses are typically modeled simply as tidal 
shear. 

How valid is this approximation? In some instances, 
higher order terms in the potential than the ~ cos 20 
term represented by external shear may be required, and 
[Kochanek (2004) describes quantitatively how to estimate 
the importance of such terms for a lens model. Addition- 
ally, the nearby environment can contribute not only tidal 
sh ear, but (projected) mass d ensity as well. As discussed 
bv lHolder & Schechteri ( l20()l . massive ellipticals are biased 
tracers of the underlying mass density; we expect massive 
ellipticals to reside in overdense environments, and indeed, 
many lens galaxies are found in poor groups or galaxy clus- 
ters. |Keeton & Zabludoff (2004), among others, point out that 
the environmental mass density, here denoted fCgrp, can lead 
to a mass sheet degeneracy in many of the quantities derived 
from strong lensing analyses. For instance, the Hubble con- 
stant inferred from time delays scales as //q °<: (1 — %p). In 
the age of precision cosmology, it is clearly important to de- 
termine the magnitude of this uncertainty in lens models. 



In principle, our theory of structure formation should al- 
ready provide an answer to this question. In the cold dark 
matter (CDM) model, the clustering of mass proceeds largely 
through the growth and merging of dark matter halos. Assess- 
ing the importance of lens environments reduces basically to 
determining, where (in which halos) do lens galaxies live? 
Placing lens galaxies chiefly in the most massive, biased ob- 
jects like galaxy clusters will of course enhance the effects of 
environment, whereas environment will be less important if 
lenses live in lower mass halos. 

We do not have complete freedom to pick the halo oc- 
cupation distribution of lens galaxies, however: there are 
now several observational constraints on models for ellip- 
tical galaxies and their hosts. First, galaxy-galaxy lensing 
(Seljak et al. 2004; Sheldon et al. 2004) provides direct con- 
straints on the average mass profiles surrounding lens galax- 
ies. Second, choosing a model for the halo occupation distri- 
bution of galaxies implicitly makes a prediction for the num- 
ber density of those galaxies, and there are now precise mea- 
surements of the number density of massive elliptical galax- 
ies (which most strong lenses are). Accordingly, we can use 
these measurements to constrain the HOD of lens galaxies, 
and thereby compute the effects of lens environments. 

The plan of this paper is as follows. In |j2] the number 
counts of Sheth et al. (2003) and the galaxy-galaxy lensing 
results of Sheldon et al. (2004) are combined, in the context 
of the halo model, to constrain the halo occupation distribu- 
tion of lens galaxies. In f|3] this HOD is used to populate an 
N-body simulation with lens galaxies, which is then exploited 
to make predictions for the effects of lens environments. 

2. HALO MODEL CALCULATION 



iHu & JainI (120031) explain in detail how the halo model may 
be used to compute the number density and average surface 
density (and shear) profiles around galaxies. As they discuss, 
the galaxy-convergence angular correlation function may be 
computed by projecting the 3-D galaxy-mass power spectrum 
using the iLimbed ( 119541) approximation. 
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Here, nv{z) is the volume number density of lens galaxies as a 
function of redshift, = / dzr^{z)-fj^nv[z) is the projected 
angular number density of lens galaxies, r{z) is the comoving 
angular diameter distance, and Dis,s the angular diameter 
distances from lens to source and from observer to source, 
respectively. 

The halo model enters by providing a prescription for the 
galaxy-mass power spectrum P^^. The ingredients are mod- 
els for the dark matter halo mass function dniJdM, the halo 
bias b{M), the halo density profile p(r), and the halo oc- 
cupation distribution of galaxies. For D M halos, we use 
the mass function and bias prescription of Sheth & Tormen| 
(0^999), the density profile of Navarro et al. ( 1997), and halo 
concentrations of Bullock et al. (2001). The halo occupation 
dist ribution (HOD) specifies how galaxies populate DM ha- 
los. |^ravteo^t^0'E004) present a model for the HOD based 
upon N-body simulations. In their model, the mean number 
of galaxies (say, brighter than luminosity L) residing in halos 
of mass M has two pieces, N{M) = Nc[M) +Ns{M), coiTe- 
sponding to central galaxies and satellite galaxies. They find 
thatA^f.(M) is well fit by a step function, A^^. (M) = ©(M-Mth) 
where the parameter Mth is the minimum mass halo capable 
of hosting a c entral galaxy with the required properties (e.g. 
luminositv). iKravtsov et alJ also find that the satellite term is 
well described by a Poisson distribution with mean 
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for M > Mth and zero otherwise. The halo mass function and 
HOD together specify the predicted galaxy volume number 
density, viz 

nviz) - J dm^N{M). (6) 

Similarly, the halo model expression for the g alaxy-mass 
power spectrum becomes, in the notation of iHu & JainI 
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y/, is the normalized Fourier transform of the DM density pro- 
file, 

yh{k,M)^M-' r"dr47zr^p{rf-^^, (11) 

and similarly y^ is the normalized Fourier transform of the 
satellite galaxy distribution within halos, here assumed to 
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Fig. 1. — Plotted are 68, 95, and 99% likelihood contours in the in, — A,, 
plane forM,h = IO'^'^/i^'Mq. 



track the mass, so that y^ = y/,. Note that the redshift- 
dependent concentration is implicitly included in the density 
profile p (r,M). The satellite mass function is given by 
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The observables to match are the number density and tan- 
gential shear profile. Eqn.|6lgives the predicted number den- 
sity, while the predicted shear profile is given by 
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Strong gravitational Isnses te nd to be massive elliptical galax- 
ies, and Sheldon et al. ^2004) report the average shear profile 
of elliptical galaxies with velocity dispersion a > 182 km/s 
observ ed in the Sloan Digital Sky Survey (SDSS). Shet h et alJ 
J2003I) have measured the velocity dispersion function of 
SDSS ellipticals, and found it to be well fit by a Schechter- 
type function. For ellipticals with a > 182 km/s, their fit 
gives a total number density of /lobs = 0.00172(Mpc//i)^^. 
Also, the mean squared velocity dispersion for this sample 
is (a^) = (216km/s)^. For a lens at z = 0.45 and source 
at Zsrc = 2, this would correspond to an Einstein radius of 
0E = 4l^{(y /c)^Dls/Ds = 0.9", in good agreement with typi- 
cal lens splittings. 

The ingredients are now in place to find the HOD param- 
eters which best fit the observed data on massive ellipticals. 
This HOD will be used in the next section to populate halos 
found in an N-body simulation, so it is important that consis- 
tent definitions of halo mass are used for the halo model cal- 
culation and the group finder, the friends-of-friends algorithm 
with linking length b — 0.2. Jenkins et al.. ( 2001) have found 
that FOF with b = 0.2 corresponds well to halo overdensity 
of 180, and that the Sheth -Tormen mass function provides a 
good fit to the mass function if the halo mass is defined not as 
the virial mass, but as M = Mi^q, i.e. the mass enclosed within 
radius rjgo for which the average interior density is 1 8 times 
the mean matter density Pm- The lBullock etai] ( l200l model 
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for halo concentrations is given in terms of the virial mass, 
not M180 which is typically ~ 1.2x larger than Myh-. Given 
the NFW density profile, it would not be difficult to translate 
from Ml 80 to Myir (Hu & Kravtsov 2003), but the dependence 
of concentration on mass is so weak (c,, °^ M^^-^^) that this 
distinction is ignored here. 

The HOD used here has three parameters: Mth, and m^. 
Figure n] shows likelihood contours in the A , —nis plane, for 
Mth = IO'^'^/z^'Mq, assuming a 10% uncertainty for the num- 
ber density nobs- Other values of Mth gave significantly poorer 
fits, so Mth = 1O'"'^/i^'M0 will be assumed throughout. As 
can be seen from the figure, a broad region in the A, — m, 
plane, centered on A, — 15, = 1.4 provides reasonable 
fits to the data. These values compare reasonably well to 
HOD parameters derived by Zheng et al. (2004) from hydro- 
dynamic simulations and semi-analytic models. The satellite 
fraction for these HOD parameters is roughly 25%, in good 
agreement with the estimates of Keeto n et al 1 ll2000l) . There 
appears to be a mild degeneracy between As and nis, in the 
sense that increasing A, (which lowers the fraction of lens 
galaxies that are satellites) may be compensated by increasing 
nis (which places satellites in more massive groups) to hold 
fixed the tangential shear profile; see fig.|2lfor an example. 
However note that large is likely disfavo red by measure- 
ments of the multiplicity function in clusters (iKochanek et al.l 

One defect which has been glossed over so far is that the 
SDSS data on elliptical galaxies used to constrain the HOD 
parameters extend to redshifts z < 0. 15, much less than typical 
lens redshifts zi ^ 0.5. In order to apply this HOD to the lens 
population, it must be assumed that there is little evolution in 
massive ellipticals between z = 0.15 and z = 0.5. 

Note that the average fCgrp may be calculated using these 
HOD parameters using just the halo model, by keeping only 
the term in Pgg corresponding to the correlation between satel- 
lite galaxies and mass in their parent halos: 

Plfik)^- [dM^^Nsy.y,. (14) 
riv J Pm uM 

This somewhat underestimates JCgip, because in this formal- 
ism halo profiles extend only out to rigo, whereas in reality 
there is mass overdensity extending well beyond this radius 
(e.g-tOuzik & Seliak 2002). With this caveat in mind, apply- 
ing ean. 1141 gives an average (fgrp « 0.02. As noted above, 
fractional errors in (e.g.) the Hubble constant will be ^((Cgrp), 
implying that estimates of the Hubble constant which ignore 
the group convergence make ^ 2% errors for most lens sys- 
tems, which is still small compared to other systematics. 

3. APPLICATION TO N-BODY SIMULATIONS 

The halo occupation distribution from the previous section 
can now be used to populate halos from an N-body simulation 
with lens galaxies, to compute the distribution of lens environ- 
ments. A small (by modern standards) dissipationless simula- 
tion with 256^ particles in a 256 /i^'Mpc box was run using 
the TPM code^ (Bode & Ostriker 2003) for a WMAP cosmol- 
ogy with Q.M ~ 0.27, Q.A = 0.73 and h = 0.72. Initial condi- 
tions were generated using the GRAFICI code^ (^ertschinger 
for a scale-invariant primordial power spectrum (n — l) 
and present-day amplitude (7g = 0.9. The particle mass for 

^ http://www.astro.princeton.edu/~bode/TPM 
' http://arcturus.mit.edu/grafic 
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Fig. 2. — The galaxy-shear correlation function for massive ellipticals. 
Data points are from Sheldon et al. 1 2004), kindly provided in electronic form 
by D. Johnston (2004, priv. comm.) The solid curves correspond to the pre- 
dicted AE = ^aitiTr) for three different halo occupation distributions, with 
As = 15, mj = 1.4; A., = 21, = 1.6; and A,, = 10.5, = 1.2. The other 
curves show the contributions to eqns. and 1101 : the galaxies' own halos 
or subhalos (dotted), group halos (dashed), and two-halo term (dash-dot). 



this simulation was m„ — 7.49 X IO'^/i^'Mq, so that halos 
of mass M ~ IO'^'^/z^'M© contained 42 particles. Since the 
simulation was dissipationless, fig = was assumed for gen- 
erating the initial conditions. 

The z = 0.45 output from this simulation was then popu- 
lated with lens galaxies using the best fitting HOD from the 
previous section, with Mth — lO'^-^/i^^M;;,, Aj ~ 13.5, and 
nis — 1-3. First, the friends-of-friends algorithm'* with linking 
length b — 0.2 was used to find halos in the simulation. Cen- 
tral galaxies were placed at the most bound particles in these 
halos. The number of satellites for each halo was determined 
by generating random numbers from a Poisson distribution 
with mean given by eqn.|5] Satellite galaxies were assigned 
the positions of randomly selected particles in each halo, to 
match the assumption above that the satellite distribution fol- 
lows the DM distribution. 

Next, convergence and shear (assuming a source redshift 
Zsrc = 2) were computed at the galaxy positions by project- 
ing the particle data along three orthogonal axes. For cen- 
tral galaxies, the contribution from the galaxies' own halos 
must be removed in order to compute solely environmental 
effects. This was accomplished by excluding, for each galaxy, 
all particles in its halo (as defined by FOF) when computing 
the projected convergence. Applying this method for shear 
would be computationally onerous, since the shear is com- 
puted by repeated Fourier transforms of the full convergence 
map. Instead, to compute the environmental shear at each 
central galaxy, comoving 100 /i^'kpc cavities were excised 
around each galaxy. This truncated particle list was then pro- 
jected to make a convergence map with comoving 31.25 /i^' 
kpc resolution, which was then Fourier transformed to pro- 

^ The implementation of FOF provided at the UW N-body shop was used, 
pittp://www-hpcc. as tro.washington.edu/tools/fof.html 
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duce shear maps of the same resolution. The resulting shear 
values at the galaxy positions were mildly dependent upon 
the precise cavity radius; using comoving 50 /i^'kpc cavi- 
ties instead of 100 /z^'kpc cavities gave shear values ^ 20% 
larger For satellite galaxies, no such excision of particles was 
required, because the satellites were placed randomly inside 
their parent groups, not on subhalos. The resulting histogram 
of fforp and 7grp values is plotted in figure |3l 

There are a few points to note. First, the mean values of 
convergence and shear are similar, with {Kg^p) ~ (l7grp|) ~ 
0.03. This agrees roughly with the simple halo model esti- 
mate for (fCgrp) given in the previous section. While the 7gip 
histogram appears well described by a lognormal distribution, 
the ffgip histogram exhibits extended tails. 

An interesting related question is that of the frequency of 
high K"gip or 72ip. In these simulations, Korp > 0.1 occurred 
with roughly ^ 6% probability, while |7grp| > 0.1 had ^ 3% 
probability. One difficulty in making these estimates is that 
the tails of the fCgip and 72ip distributions are likely overrepre- 
sented, due to projections. Galaxies which project near unre- 
lated massive groups will appear to have higher fCgrp and 7grp 
than is truly representative for their own environments. Sim- 
ilarly, galaxies which project near voids will have overly low 
fCgrp. Of course, from the viewpoint of estimating errors on 
derived quantities like the Hubble constant, such projections 
are just as important to consider as the immediate environ- 
ment of the lenses. 

One last point to consider is that selection effects can mod- 
ify the probability distributions for JCgrp and 7grp for actual 
lenses compared to the distributions computed here. For ex- 
ample, the strong lensing cross section scales as (1 — ffgip)^^, 
so in principle the histograms of convergence and shea r 
should be weighted by this factor (iHolder & Schechteil2003l) . 
Reweighting the distributions (excluding the tiny fraction of 
galaxies with Kgrp > 0.5), the average (fCgrp) increases to 
0.037, while (7gip} increases to 0.033. Another possible se- 
lection effect is that ellipticals with higher velocity dispersion 
(and hence larger lensing cross section °^ 9^°^ cr^) will tend 
to live in more massive, more biased halos, and so the K and 
7 histograms should be weighted by this factor as well for 
lenses. Unfortunately, this requires a more detailed halo oc- 
cupation distribution than has been used here. 

4. DISCUSSION 

This paper has presented a calculation of the effects of the 
environments surrounding lens galaxies. Lens galaxies were 
placed in an N-body simulation using a halo occupation dis- 
tribution calibrated to match both the number density and tan- 
gential shear profiles observed for massive elliptical galax- 
ies in the Sloan Digital Sky Survey. The SDSS results are 
measured for galaxies at low redshift, z < 0.2, while lens 
galaxies are typically at higher redshift, z ~ 0.5, so the re- 
sults presented here will be meaningful only if there is little 
evolution in elliptical galaxies between z « 0. 15 and z ~ 0.45. 
The expected average values for convergence and shear are 
(jfgrp) ~ (|7grp|) ~ 0.03, with a spread of ^ 0.6 and 0.35 
in logjQ for K and 7 respectively. Values much higher than 
the mean appear relatively rarely, with fCgrp > 0. 1 occurring 
with ~ 6% probability and 7grp > 0.1 occurring with 3% 
probability. 

Previous work has considered related topics . Both 
iKeeton eta f ( 1997) and Holder & S chechter (2003. HS) have 
estimated the magnitude of tidal shear produced by mate- 
rial correlated with lens galaxies. This paper has followed 



an approach quite similar to that developed b y HS. The re- 
sults presented here appear similar to those of iKeeton et all 
however HS find somewhat higher environmental shear. This 
appears to be due to a different HOD; HS use massive el- 
liptical galaxies produced using the semi-analytic recipes of 
■Kauffmann e t al. ( 1999). As they note, their lens galaxies tend 
to live in very massive halos, with roughly ~ 30% residing 
in clusters more massive than IO^^Ii^^Mq. Accordingly, the 
number density of massive elliptical galaxies predicted by HS 
is lower than the values computed here by roughly an order 
of magnitude. Further observations of the environments sur- 
rounding lens galaxies (e.g. |Fassnacht et al. 2004) can help 
determine whether lenses typically do reside in such massive 
halos. 

K.\\'ii:l\. Keeton & Zabludoff ('2004') have investigated the 
importance of the group convergence on strong lenses. They 
construct a mock group including several supercritical galax- 
ies separated by a few Einstein radii. By explicit construc- 
tion, they showed that in cases like these, ffgrp can approach 

0.15 — 0.2, which would dramatically skew the results of 
analyses ignoring group convergence. Clearly, for systems 
like these, it is necessary to model such close neighbors us- 
ing both convergence and shear, and indeed in systems like 
PG 1115H-080 the nearby groups are modeled with isother- 
mal spheres or NFW profiles (e.g. iKochanek & Dalal"2004t 
Treu & Koopmans 2002). For the majority of systems, how- 
ever, the environment is treated merely with tidal shear, and 
the results of this paper indicate that such a treatment is usu- 
ally adequate. If typical (Cgrp at lens galaxies is ^ 3%, then 
errors in the Hubble constant will similarly be at the ^3% 
level, which is still unimportant compared to the > 10% sys- 
tematic errors typical of lensing Hq's. 

One remaining puzzle from this work is the origin of the 
large shears measured at lens galaxies. As discussed by 
'Keeto n et alJ ( Il997h . successful models of strong lenses typi- 
cally require external shear at the level of 7^ 10— 15%, con- 
siderably higher than the 7grp « 3% found here. Such shear 
cannot be explained by uncorrelated structure (e.g. line-of- 
sight projections), which should produce shear of order the 
cosmic shear, ^ 1 — 2% on these scales. HS have argued that 
such strong shear naturally occurs for lenses living in mas- 
sive systems like gal axy clusters w here JCv = Psfs/^ah « 0. 15. 
On the other hand, Xeeton et alJ have suggested that much 
of this so-called external shear originates in the halos of the 
lens galaxies themselves. To support this idea, they note that 
the total quadrupole at the images aligns quite well with the 
galaxy orientation, while misalignments larger than observed 
would arise if 10 — 15% external shear were completely un- 
correlated with the galaxy orientation (see^ochanek 2004 
for an updated version of this). As noted by Holder (2003, 
priv. comm.), however, this argument has difficulty. Halo ma- 
terial stratified on ellipses with the same ellipticity as the mass 
interior to the Einstein ring will produce no effect on the lens 
images, as shown by Newton. For the halo to produce this 
shear, it must either have principal axes substantially twisted 
relative to the galaxy principal axes, or it must have a very 
different ellipticity. To illustrate, imagine taking an isother- 
mal ellipsoid with Einstein radius rg, and beyond some radius 
ro, either twisting the ellipsoidal axes by some angle AO or 
changing the ellipticity by Ae. Then the change in shear at 
the Einstein radius, in the limit ro ^ ^Ej will be of order 

A7- (Ae + 2esin2A0), (15) 
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Fig. 3. — Histograms of (left) (fgip and (right) |7grp| for the z = 0.45 simulation, assuming Zsrc = 2. The HOD used here had M,i, = 10'^ '/i-'Mg, A., = 13.5, 
and mj = 1.3. 



where typical lens Einstein r adii are ~ 2reff and typical el- 
lipticities are e 0.3. Hoekst ra et alJ (12004) have found that 
galaxies and their extended halos align quite well, with halo 
ellipticities on average > 0.7 x the observed galaxy ellipticity, 
which limits A7< Q3er^/2nrQ. This casts doubt on a halo 
origin for observed external shear; ro ^ rg would be required 
to produce enough shear in this way. Therefore, dramatic 
changes in the density structure must occur in lens galaxies 
just outside their Einstein radii, near the transition between 
baryon-domination to DM-domination, for the lenses them- 
selves to account for the observed shear. On the other hand, 
halo substructure conceivably could account for some of the 



unexplained shear (ICohn & Kochanekl2004) . 
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